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Atmospheric C0 2 concentrations appear to have been considerably higher than modern levels during much 
of the Phanerozoic and it has hence been proposed that surface temperatures were also higher. Some studies 
have, however, suggested that Earth's temperature (estimated from the isotopic composition of fossil shells) 
may have been independent of variations in atmospheric C0 2 (e.g. in the Jurassic and Cretaceous). If large 
changes in atmospheric C0 2 did not produce the expected climate responses in the past, predictions of 
future climate and the case for reducing current fossil-fuel emissions are potentially undermined. Here we 
evaluate the dataset upon which the Jurassic and Cretaceous assertions are based and present new 
temperature data, derived from the isotopic composition of fossil brachiopods. Our results are consistent 
with a warm climate mode for the Jurassic and Cretaceous and hence support the view that changes in 
atmospheric C0 2 concentrations are linked with changes in global temperatures. 

During the Phanerozoic, climate has alternated between states characterised by widespread glaciation (ice- 
houses) and largely ice-free hothouse (or warmhouse) conditions 16 . Quasi-periodic icehouse-hothouse 
cycles of varying length have been proposed 79 and many 2,8,10 have suggested that that these major changes 
in Phanerozoic climate are linked to variations in atmospheric C0 2 . A number of studies have, however, cast doubt 
on the C0 2 -climate link 4 and suggest instead that the galactic cosmic ray flux is the main driver of Phanerozoic 
climate change 511 . Although not fully accepted by many 3,12,13 , if correct, then this view has significant implications 
for our predictions of future climate change and the case for reducing anthropogenic carbon emissions 6 . 

The basis for this assertion is a putative temperature record 51415 based upon the oxygen isotopic composition of 
well-preserved marine fossils (principally brachiopods, belemnites and planktonic foraminifera). This temper- 
ature record 5 shows two mismatches between temperature and atmospheric C0 2 : during the Tate Ordovician 
glaciation (—445 million years ago) and the Early Jurassic to Early Cretaceous interval (—180-100 million years 
ago), both of which coincided with model 216 and proxy indications 10,18 of high levels of atmospheric C0 2 . The Late 
Ordovician paradox has been particularly well studied 17-20 . A lower solar luminosity at this time suggests that the 
C0 2 threshold for initiating a glaciation may have been higher 18 . Moreover, if the cooling pulse that led to ice 
volumes to expand to a similar size as those of the Pleistocene was relatively short-lived 19 , atmospheric C0 2 
proxies may not have recorded this brief reduction in C0 2 levels. Furthermore, it has also been suggested that 
tropical ocean temperatures in excess of 30°C during the Late Ordovician coexisted with substantial south polar 
ice sheets 19 . Significant differences in the geosphere, e.g. the distribution of the continents, and biosphere, e.g. the 
absence of significant terrestrial vegetation and calcareous plankton, during the Late Ordovician suggest that the 
links between C0 2 and climate may well have operated differently to more recent times. 

The Early Jurassic to Early Cretaceous anomaly presents a more critical test of the relationship between C0 2 
and temperature change, given that Earth's palaeogeography and biota were more similar to those of the present- 
day. Although the Mesozoic is not a direct analogue for future greenhouse warming, such warm intervals in Earth 
history provide important insights into processes operating in the climate system 21 . Isotope derived temperature 
curves for this interval 4 suggest significant cooling at this time, whereas C0 2 proxy and model data typically 
indicate that atmospheric C0 2 reached peaks not recorded since the Devonian 2 . Here we present extensive 
temperature data, derived from the isotopic composition of fossil brachiopods in order to evaluate the hypothe- 
sized Jurassic-Cretaceous cooling trends and re-assess Phanerozoic trends. 

Results 

Over 170 new oxygen isotope measurements from Mesozoic calcitic brachiopods are presented, combined with a 
further 300 brachiopod-derived measurements from published data 22 25 (Figs. 1, 2) in order to test the robustness 
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Figure 1 | Evolution of Mesozoic-Cenozoic 5 ls O values. Smoothed curves (red curve 12,13 and brachiopod-foraminifera only isotope data, blue curve) are 
generated using Kernel regressions with a bandwidth (h) of 1 0 Myr. Abbreviation: Trias = Triassic. Grey dots = brachiopod data, red circles = belemnite 
data and green circles = planktonic foraminifera 12,13 and black circles = brachiopod data compiled in this study. 



of cooling events suggested in previous studies 4,5 , and test whether 
galactic cosmic ray flux is potentially a driver of Phanerozoic climate 
change. The new dataset (see SI) covers the Jurassic to Cretaceous 
time interval. Almost all shells are derived from shallow epicontin- 
ental seas of the European region representing the tropical to sub- 
tropical Tethys. The data are constrained by ammonite biozone 
resolution (—0.4 to 1 Myr on average) allowing the generation of a 
reliable temporal series calibrated to numerical ages 26 . Our cathodo- 
luminescence (CL) and petrographic analyses (SI) of the calcitic 
shells demonstrate excellent ultrastructure preservation of the ana- 
lysed material. The trace element concentrations of each of the bra- 
chiopod shells corresponds to the results of the petrographic and CL 
analyses, whereby relatively low concentrations of Mn and Fe also 
argues for well-preserved fossils 3,14-15 . 

Discussion 

Our new brachiopod-derived 8 ls O isotope data from the Jurassic and 
Early Cretaceous are consistently more negative than coeval pub- 
lished data 1415 (Fig. 1). The most likely explanation is that the pub- 
lished data derive from different fossil taxa; the Jurassic and Early 
Cretaceous portion of the published database is dominated by data 
derived from belemnites (Cephalopoda, Coleoidea). Studies 24,25-27 
have also shown that belemnites typically display relatively more 
positive 8 ls O values than brachiopods; a difference possibly related 
to a combination of factors. Firstly, belemnites are interpreted as 
migratory nektobenthic organisms that were hence able to tolerate 
deeper and cooler waters, just as extant coleoids do 28 . Secondly, they 
may not have precipitated their calcite in isotopic equilibrium with 
seawater 25 . Furthermore, belemnites within the database are derived 
from a range of low, mid and high palaeolatitudes 15 whereas the 
brachiopod data are derived predominantly from low (tropical and 
subtropical) palaeolatitudes 14,15 . 

In contrast, the Late Cretaceous-Cenozoic proportion of the pub- 
lished database 14,15 consists of isotope measurements derived largely 
from planktonic foraminifera. Where our new data and published 
brachiopod data 24 are coeval with these planktonic data they show a 
good correspondence. The shells of planktonic foraminifera and 
brachiopods of similar age record comparable 5 ls O signatures, 



whereas coeval belemnites do not. Has, therefore, the inclusion of 
belemnite data in previous datasets 14 resulted in erroneous temper- 
ature estimates for the Early Jurassic-Early Cretaceous interval? 

As we have excluded a diagenetic influence (and assuming shell 
calcite is precipitated in isotopic equilibrium with seawater 22 25 ), the 
8 18 0 composition of the brachiopods and planktonic foraminifera 
must reflect both the isotopic composition of the seawater in which 
the organisms lived and the water temperature at the time of shell 
growth. The isotopic composition of the seawater reflects changes in 
the size of ice-sheets, (ice-sheet growth preferentially removes 1<s O 
from seawater), the pH of seawater and possibly a long term change 
in seawater caused by tectonic processes 29 . Recent research using the 
'carbonate clumped isotope' method 3,19 is, however, inconsistent 
with such a secular variation in the 8 ls O of sea water. For this reason 
we have not adjusted the isotope data to accommodate for this 
(Figs. 1, 2), and the curve is interpreted in terms of temperature 
change. Given the potential problems with including belemnite data, 
we have combined our new data with the brachiopod- and forami- 
nifera-derived data 14,15 only, without including any belemnite data, to 
produce a novel, high resolution Mesozoic temperature curve. 

When the critical Early Jurassic to Early Cretaceous interval is 
examined, it is clear that the new curve of mean temperature diverges 
significantly from the previous one (Fig. 1). The maximum difference 
in S 18 0 is 2.5 %o, equivalent to ca. 8°C mean temperature change. We 
suggest the inclusion of belemnite data 4 skews interpretations by 
creating the impression that temperatures were some 8°C cooler than 
the new data suggest. 

The excellent agreement between our new Mesozoic brachiopod 
data and the published data derived from coeval planktonic forami- 
nifera enables us to confidently link the published Palaeozoic dataset 4 
(derived almost entirely from brachiopods) with Cenozoic data 
(derived almost entirely from planktonic foraminifera), providing 
a new Phanerozoic 8 ls O (temperature) curve that is free of biases 
related to the use of belemnite data in the Mesozoic. To test the 
relationship between temperature and C0 2 at various temporal 
scales through the Phanerozoic, we use the "SiZer" statistical 
approach (Significance of Zero Crossings of the Derivative 30 ). This 
method is based on the construction of curves fitting time series 
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Figure 2 | Evolution of Phanerozoic 5 18 0 values. (A): Smoothed brachiopod-foraminifera only curve generated using Kernel regressions with a 
bandwidth (h) of 10 Myr. Abbreviations: Camb = Cambrian; Ord = Ordovician; Sil = Silurian; Carbon = Carboniferous; Perm = Permian; Trias = 
Triassic. Green shaded boxes represent glaciations and icehouses 12 . (B): SiZer map of the isotope data derived from brachiopods and foraminifera only. 
Blue shading indicates a statistically significant increase in isotopic values (i.e. cooling), and the confidence interval is completely above a slope of 0. Red 
shading indicates that the confidence interval is completely below a slope of 0, and there is a significant decrease (i.e. warming). Purple shading indicates 
where the confidence interval for the derivative contains 0. White fill indicates where the data are too sparse. The lower panel C shows the atmospheric C0 2 
concentration record as determined from multiple proxy reconstructions 16 (green curve), and C0 2 estimates (yellow curve) from GEOCARB IIP. 
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using different levels of temporal smoothing (h). The first derivatives 
of each curve (i.e., the slopes) are computed with their 95% confid- 
ence intervals, allowing the signs of derivative estimates to be stat- 
istically tested. The results of these multiple tests reported in the 
form of SiZer maps (Fig. 1) enable significant features to be identified 
at different time scales. SiZer has been successfully applied to 
studies of the Quaternary 31 and of the Jurassic 27 , but this is the first 
Phanerozoic-scale analysis. 

The long-term, first-order Phanerozoic trend is one of overall 
cooling (i.e. increase in 8 ls O) from the Cambrian/Ordovician bound- 
ary to the present day, which matches the first-order decline in C0 2 
over the same interval (Fig. 2). Our SiZer analysis of the Palaeozoic 
record indicates two major cooling episodes (through the Ordovician 
and from the mid-Devonian to end-Carboniferous), both of which 
culminate with well-known glacial maxima 9 . Significant repeated 
shorter-term trends are apparent within the Ordovician and 
Silurian, a time characterised by warming and cooling/glacial and 
interglacial cycles 19,20 . We are therefore confident that our approach 
will provide a robust means to evaluate Mesozoic-Cenozoic trends. 

Our SiZer analysis reveals warming during the Permian corres- 
ponding closely to rising C0 2 levels, reaching a peak around the 
Permian/Triassic boundary 32 . Following the Permian/Triassic tem- 
perature maximum, our new Mesozoic brachiopod S ls O record sug- 
gests an overall cooling trend through most of the Triassic, mirroring 
the C0 2 estimates from the GEOCARB III model 2 and the 
GEOCLIM 32 and COPSE models 33 . There is less agreement between 
our data and other proxy C0 2 estimates 18 in particular those derived 
from paleosol carbonates 10 . Through the Jurassic, the long term tem- 
perature trend is one of relative stasis. Evidence for long lived and 
extensive ice sheets during this time is lacking, although short tran- 
sient events may characterise this time 8,34 . The SiZer analysis only 
records one cooling episode in the Jurassic. Indeed, for this interval a 
wealth of data is available suggesting that high latitude warmth is the 
norm 1,8,18 . The apparent long term climatic stasis continues through 
to the Late Cretaceous (ca. 95 Ma), when a return to increasing 8 ls O 
values indicate the onset of another cooling trend, consistent with 
temperature predictions from the COPSE model 33 . The transition 
from equable temperatures to long-term cooling is observed as 
C0 2 levels fell below 1000 ppm. The long-term trend of climatic 
stability for at least 100 million years (largely encompassing the 
Jurassic and Early Cretaceous) recorded in our brachiopod 8 ls O data 
is also reconcilable with the C0 2 predictions for this interval, 
although there is some disagreement between model and proxy esti- 
mates of atmospheric C0 2 levels 10,16,18 (Fig. 2C). 

Our new brachiopod-derived 8 ls O isotope curve overturns the 
earlier and contentious suggestions that during the Early Jurassic 
to Early Cretaceous Earth experienced a cool climate mode. 
Inclusion of belemnite 8 18 0 data is shown to bias the temperature 
interpretation towards cooler estimates. Whilst the data are consist- 
ent with two icehouse-hothouse intervals lasting tens of millions of 
years within the Palaeozoic, the Mesozoic-Cenozoic shows a differ- 
ent pattern. The lack of a prolonged cooling episode leading to glacial 
conditions may be a consequence of elevated atmospheric C0 2 
coupled with changing continental positions and higher solar lumin- 
osity. The Permo-Triassic hence appears as a major turning point in 
the climatic, as well as biotic 35 , history of the Earth. Carbon cycle 
modelling estimates of C0 2 changes over the Mesozoic, in particular, 
correlate well with these changes in temperature. Together our obser- 
vations support the assertion that atmospheric C0 2 is consistent with 
patterns of past global climate change on a Phanerozoic time scale, 
and serves to invalidate claims against a long-term C0 2 -climate link. 

Methods 

For stable isotope ratios we use the conventional terminology and 5-notation: 5 18 0 — 
(Rsample/Rstandard-1) * 1000. Rsample and Rstandard are the ls O/ 16 0 ratio in the 
sample and standard, respectively. Positive or negative 8 !s O values mean that the 
sample has more or less la O in the standard in parts per thousand (permil, %o). The 



PDB standard, is a Cretaceous marine fossil (Belemnitella americana), and has now 
been exhausted and replaced by the V-PDB standard. Stable isotopes were deter- 
mined on a VG Instruments Optima Isotope Ratio Mass Spectrometer (at the 
University of Plymouth) using 200-300 ug carbonate. Isotopic results were calibrated 
against NBS-19. Reproducibility for both S l8 0 and 5 13 C was better than 0.1%o, based 
upon multiple sample analysis. Trace element geochemistry in conjunction with 
optical petrography and cathodoluminescence (CL) was used to determine the state 
of preservation of each of the brachiopods analysed. Cold cathode CL analysis was 
undertaken using a CITL CL MK 3A, mounted on a Nikon petrological microscope. 
Sub-samples for chemical analysis weighing 1-2 mg were dissolved in nitric acid and 
analysed using a Varian 725-ES ICP-AES. Based upon analysis of duplicate samples, 
reproducibility was better than +/— 3% of the measured concentration of each 
element. 
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